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ABSTRACT 

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes 
O , of galactic gravitational potentials since two rotation curves can be measured in two 

I perpendicular planes. Observational studies have uncovered systematically larger ro- 

O tation velocities in the extended polar rings than in the associated host galaxies. In the 

■4-^ dark matter context, this can only be explained through dark halos that are system- 

^ atically flattened along the polar rings. Here, we point out that these objects can also 

I I be used as very effective tests of gravity theories, such as those based on Milgromian 

dynamics (MOND). We run a set of polar ring models using both Milgromian and 
^~^ Newtonian dynamics to predict the expected shapes of the rotation curves in both 

^ planes, varying the total mass of the system, the mass of the ring with respect to the 

^~~' host, as well as the size of the hole at the center of the ring. We iind that Milgromian 

^2J dynamics not only naturally leads to rotation velocities being typically higher in the 

~^ l extended polar rings than in the hosts, as would be the case in Newtonian dynamics 

without dark matter, but that it also gets the shape and amplitude of velocities cor- 
\l rect. Milgromian dynamics thus adequately explains this particular property of polar 

^--^ ring galaxies. 
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1 INTRODUCTION at which the mass discrepancy between baryonic and dy- 

namical mass appears in the standard picture, the tran- 
Assuming General Relativity to be the correct description ^.^.^^ ^^^^^^^ acceleration between dark-matter dominated 

of gravity at all scales, data ranging from the largest scales , , , • i j i • / -iu- at i • 

, . . \ ■ ^^" baryons-dommated galaxies (witnm Newtonian grav- 



^ (e.g., the Cosmic Microwave Background) to galactic scales .^^-^^ ^^^ .^ ^^g^^^ ^ ^^.^.^^^ ^^^^ ^^^^^^^ ^^^^^^ ^^^ ^.^^ 

C^ can be interpreted as a Universe dominated by dark energy ^^^^.^.^^ (Famaey & McGaugh 2012) . These independent oc- 

and dark matter. The nature of these is among the most r iiii jij-iu^jj 

° curences oi ao are not at all understood m the standard con- 

challenging problems of modern physics. While dark energy ^^^^^ whereas, surprisingly, all these relations can be sum- 
is generally assumed to be a non-vanishing vacuum energy ^^^.^^^ ^^ ^^^ empirical formula of Milgrom (1983). For 
represented by a cosmological constant A in Einstein's equa- ^^.^ ^^^^^^^ ^^ g^ ^^1^^^ ^^^^^.^^ ^^^^^^^ ^^^ ^^^^^ q^^^^^ 
tions, the currently most favored dark matter candidates are ^^j^^ ^^ ^^^ ^^^ ^^^^ ^^1^ j^^^^^^^ q 9 ^^^ ^5 ^ ^q-io^ ^-2^ 
neutral fermionic particles, which condensed from the ther- ^^^ ^^^^ ^ ^^^^^ .^ ^g^^^^^ ^^ ^^^^^^ ^^^^^.^^ ^^^^^ ^^^^ 
mal bath of the early Universe (Bertone et al. 2005; Strigari ^^ g^^^^ ^.^^ ^ ^.^^^^ ^^1^^ (^^^^g^ ^^ ^1 2011). We choose 
2012), known as "cold dark matter" (CDM) particles. g^^^ ^g^ ^^^.^^ ^^1^^ ^^ ^ ^2 x lO-^Vs-^, as per Gentile 

On galaxy scales, predictions of this concordance cos- ^ i foOlbl 
mological model (ACDM) are difficult to reconcile with 

observations (Disney et al. 2008; Peebles & Nusser 2010; The success of Milgrom's empirical formula lends weight 

Kroupa et al. 2010; Kroupa 2012; Kroupa et al. 2012). For to the idea that the gravitational field in galaxies can be 

instance, many observed scaling relations (see Famaey & described by Milgromian Dynamics (also known as Modi- 

McGaugh 2012, for a review) involve the universal appear- fied Newtonian Dynamics or MOND). Milgromian dynamics 

ance of an acceleration constant ao ~ A^" ~ 10^^"ms^^ ~ naturally explains the intimate relation between the distri- 

3.6 pc Myr~'^, whose origin is unknown in the standard con- bution of baryons and the gravitational field in galaxies, and 

text. For instance, this constant defines the zero-point of explains all the aforementioned occurrences of oo in galactic 

the TuUy-Fisher relation, the transition of the acceleration dynamics without any fine-tuning. Given the predictive na- 
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ture of Milgromian dynamics on galaxy scales, it is of great 
interest to test whether the formula can explain all probes 
of galactic gravitational potentials, beyond spherical and ax- 
isymmetric systems where it has mostly been tested up to 
now. 

Polar ring galaxies (PRGs) are non-axisymmetric sys- 
tems featuring an outer ring of stars and gas rotating over 
the poles. The host galaxy is usually characterized by a com- 
pact bulge and a small bright gas-poor disk, while the gas- 
rich polar structure has photometric properties roughly sim- 
ilar to those of gas-rich spirals (e.g., Whitmore et al. 1990). 
The observer can typically measure two perpendicular rota- 
tion curves (Schweizer et al. 1983; Sackett & Sparke 1990; 
Sackett et al. 1994; Reshetnikov & Combes 1994; Combes 
& Arnaboldi 1996; lodice et al. 2003, 2006; lodice 2010), 
one in the host, often by deriving an asymmetric-drift cor- 
rected rotation curve from the observed stellar kinematics 
(see, e.g., Combes & Arnaboldi 1996), and one in the po- 
lar ring, by directly measuring the velocity of the HI gas. 
This makes PRGs ideal test objects for gravity theories, be- 
cause any given theory of gravity then has to explain two 
rotation curves in two perpendicular planes, both derived 
from the same baryonic mass density distribution. Inter- 
estingly, observational studies (lodice et al. 2003; Moiseev 
et al. 2011) consistently show rotational velocities in the po- 
lar rings to be systematically larger than in the hosts. These 
observations may only be explained in the standard context 
by dark halos systematically flattened along the polar rings 
(see lodice et al. 2003). In any case, given these specific 
observational properties of PRGs, it is of great interest to 
investigate whether the general predictions of Milgromian 
dynamics for such objects would conform with these obser- 
vational properties, namely whether larger velocities in the 
extended polar rings than in the hosts are a generic predic- 
tion of Milgromian dynamics, by exploring a wide range of 
baryonic mass distributions. 

In Sect. 2, we recall the basics of Milgromian dynam- 
ics and the specific quasi-linear formulation we are dealing 
with. We then present a grid-based prescription to solve the 
modified Poisson equation (Sect. 3) and an iterative method 
to find rotation curves of non-circular orbits (Sect. 4), and 
apply it to a set of models in Sect. 5. Results are presented 
and discussed in Sect. 6 and we conclude in Sect. 7. 



2 MILGROMIAN DYNAMICS 

In recent years, a plethora of generally covariant modified 
gravity theories have been developed, yielding a Milgromian 
behavior in the weak-field limit (Famaey & McGaugh 2012). 
One such recent formulation (Milgrom 2009) has a non- 
relativistic quasi-static weak field limit, for a specific given 
set of parameters, yielding the following Poisson equation: 



V^$ = AnGpb + \7-[i' (! V<^i/ao) V0] , 



(1) 



where $ is the total (Milgromian) potential, pb is the bary- 
onic density, 4> the Newtonian potential such that \7^(j> — 
AnGph, and where i^(a;) — >■ for x ^ 1 and t^ix) — > x^^" 
for X ^ 1. One family of functions that fulfills the definition 
of !^(x) (see, e.g. Famaey & McGaugh 2012) is 



Hereafter, when not stated otherwise, we use n = 1, a func- 
tion which is known to reproduce well the rotation curves of 
most spiral galaxies^ (Gentile et al. 2011). 

This means that the total gravitational potential $ — 
(jy + ^ph can be divided into a classical (Newtonian) part, 0, 
and a Milgromian part, $ph. The matter density distribu- 
tion pph that would, in Newtonian gravity, yield the addi- 
tional potential $ph, and therefore obeys V^"l?ph = 47rGpph, 
is known in the Milgromian context as the "phantom dark 
matter" (PDM)^ density. 



Pph 



V-[//(|V(^i/ao)V0] 
4tvG 



(3) 



This is the density of dark matter that would boost the New- 
tonian gravitational field to give precisely the same effect 
as the boost of gravity predicted by Milgromian dynamics. 
For a disk galaxy, it will typically resemble a round isother- 
mal halo at large radii, but exhibits an additional disk of 
phantom dark matter aligned with the baryonic disk (but 
with a different scale length and height) most prominent at 
smaller radii (Milgrom 2001). At each spatial point, pph is a 
non-linear function of the Newtonian potential. As the non- 
linearity of Eq. 1 is only present on the right-hand side of 
the equation, it is called the quasi-linear version of Modi- 
fied Newtonian dynamics (Milgrom 2010) whereas in older 
versions of Milgromian dynamics theories the Laplacian op- 
erator on the left-hand side was replaced by a non-linear 
one (Bekenstein & Milgrom 1984). In the next section, we 
present a grid-based prescription to calculate the PDM den- 
sity. We will then be able to compute the rotation curves 
from the velocity of the closed orbits crossing the planes of 
symmetry of the non-axisymmetric system (i.e. the plane of 
the host galaxy and of the polar ring). 



u(x)= [l + (l + 4x-")^/V2] 



l/n 



1. 



(2) 



3 GRID-BASED CALCULATION OF THE 
PHANTOM DARK MATTER DENSITY 

The PDM density that would source the Milgromian force 
field in Newtonian gravity is defined by Eq. 3 and can be 
calculated from the known classical (Newtonian) potential 
0. To evaluate this term, we devise a numerical, grid-based 
scheme that calculates Pph from any (discrete) Newtonian 
potential (j)^'^'^ (see Famaey & McGaugh 2012, Eq. 35; Angus 
& Diaferio 2011; Angus et al. 2012). 



^ Nowadays, galaxy data still allow some, but not much, wiggle 
room on choosing the interpolating function u{x) (Eq. 2): they 
tend to favor the n = 1 function from the family used here; some 
interpolation between n = 1 and n = 2; or functions from other 
families which actually reduce, for accelerations typical of galax- 
ies, to the n = 1 case used here. See, e.g.. Sect. 6.2. of Famaey & 
McGaugh (2012) for a review. 

^ In the Milgromian context, PDM is not real matter but a nu- 
merical ansatz which helps to compute the additional gravity pre- 
dicted by Milgromian dynamics and gives it an analog in Newto- 
nian dynamics. 



The discrete form of Eq. 3 then reads 

ttG h^^ V v y Bx 



Pph 



4tvG I 



(4) 
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with h being the constant one-dimensional grid step size. 
The whole discretization scheme is illustrated in Fig. 1. Note 
that this equation was first derived for older versions of Mil- 
gromian dynamics by Brada & Milgrom (1999) and Tiret & 
Combes (2007) and we derive it for the quasi-linear formu- 
lation here (see also Angus & Diaferio 2011; Angus et al. 
2012). The function iy{x) is evaluated at the points marked 
by squares in Fig. 1. The gradient of in «/ (| V(jf)|/ao) at the 
point Bx has been approximated by 



V(/>: 




(5) 



Having evaluated this field, the Newtonian Poisson equation, 
V^$ = AttG {ph + Pph), can be solved to find the effective 
Milgromian force field. This can be done using the same 
grid. In the present paper, grids with a resolution of h — 
0.23 kpc are used to calculate the rotation curves (see next 
section) from 2 to 16 kpc, and grids with h — 0.47 kpc for 
radii larger than 16 kpc. The resolution was chosen such that 
it is sufficiently fine that the form of the rotation curves does 
not change if the resolution is further increased. 



4 CALCULATION OF ROTATION CURVES 

From a given Milgromian potential <1> (see next section for 
detailed description of PRG models), rotation curves (or 
rather, their non-circular equivalent in non-axisymmetric 
configurations) can be calculated. In an axisymmetric po- 
tential, the circular rotation velocity v(r), which results 
in closed orbits with radius r, readily follows as v{r) — 
\/—r ■ d$/dr in the plane of the galactic disk. This equa- 
tion however loses validity in a non-axisymmetric potential 
like the one of a PRG, because the closed orbits are gen- 
erally not circular. The existence of two massive systems 
in perpendicular orientations means that circular orbits do 
not exist in either system, neither the equatorial nor the po- 
lar one: in each plane, the potential well corresponding to 
the other perpendicular system produces the equivalent of 
a (non-rotating) bar along the line of nodes. In that case, it 
becomes necessary to obtain the velocities in the disk and 
polar ring in a more general way. In the present work, an 
iterative method is applied: test stars are shot through the 
galactic potential, which is computed numerically from an- 
alytical density distributions following the prescription in 
Sect. 3. The initial velocity (perpendicular to the radius) of 
these test particles is adjusted until a closed orbit is found. 
The orbit is integrated using the simple leapfrog integration 
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Figure 1. Illustration of the discretisation scheme in the x-y 
plane, referring to Eq. 4. The grid has A^'' nodes (i,j,k) with 
i,j,k G {1 . . . N} that are separated by a constant grid step size 
h. The values of !^(x) are evaluated at the points Ajj and B^. in 
fc-direction using Eq. 5. 



Line of sight 




Figure 2. Closed orbits (black, solid lines) within the potential 
of the benchmark model (Seq. 1, Mpj^ = 1.45 M^i^k) in the plane 
of the host galaxy. The polar ring, which is located in the y-z 
plane, is illustrated by the thick, grey line. The blue, dashed line 
is a circle which, by comparison, demonstrates the non-circularity 
of the closed-loop orbits. The major axis of the eccentric orbits 
points in the x-direction, because the test particles orbiting in 
the galactic disk (x-y) plane "fall" through the polar ring (y-z) 
plane, i.e. they feel a stronger acceleration in x than in y direction. 
However, to fulfill closed orbits, the oscillation period in both 
directions must be the same, which means that the oscillation 
amplitude in the x direction must be larger (major axis) than 
that in the y direction (minor axis). The rotation velocity thus 
is minimal at point A (along the line of sight) and maximal at 
point B. 



scheme. Typical closed orbits in a PRG potential are shown 
in Fig. 2 (detailed description of the model in the next sec- 
tion). 
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Figure 3. The polar ring density profile is made from the dif- 
ference of two Miyanioto-Nagai density profiles (pi and p2). Its 
density is zero at the centre. This example is a demonstration 
of a ring with a total mass M = 9.5 X 10^ Mq, a scale height 
hz = 0.3 kpc and scale radii hri = 6.8 kpc and h^^ = 5.95 kpc 



5 MODELS 

In order to explore the consequences of Milgromian dynam- 
ics for the rotation curves in polar rings, we start from 
a benchmark model adopted from Combes & Arnaboldi 
(1996), which represents a prototypical example of polar ring 
galaxy (NGC 4650 A). From this model, we will construct a 
Milgromian potential in which the orbits of test particles 
will be computed. The host galaxy is made of a small Plum- 
mer bulge (Plummer 1911) weighing My, = 0.2 x IO'^Mq, 



with a Plummer radius r^ = 0.17 kpc, 



-«^(S|)"+ 



-5/2 



(6) 



and of a Miyamoto-Nagai disk (Miyamoto & Nagai 1975) 
with disk mass Md = 11 x 10^ Mq, scale-length hr = 
0.748 kpc and scale-height h^ = 0.3 kpc, 

Pd(i?,^)=(^)x 



(7) 



To this parent galaxy a polar ring of stars and another 
one of gas is added^. Each ring is built by the differ- 
ence of two Miyamoto-Nagai density distributions of scale 
height /iz = 0.3 kpc, and with scale radii h,.-^ and h,-^ (see 
Fig. 3 for an illustration, where one sees that the bary- 
onic density is precisely zero at the center and positive 
elsewhere). The masses of these two disks are chosen such 
that their difference equals the total mass of the ring and 
that the central mass density of the ring is zero. The stel- 
lar ring weighs 9.5 x 10^ Mq, and has ft^* = 6.8 kpc and 
/irj = 5.95 kpc, while the gaseous ring weighs 6.4 x 10® M©, 
and has ftff'' = 15.3 kpc and /if^'' = 3.4 kpc. The parame- 



^ We do not consider the possibility of two gaseous disks to avoid 
orbits crossing, unless there is a very small gas disk and a much 
larger gaseous polar ring which never intersect. In our models, we 
make no distinction between gas and stars in the host disk. 
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Figure 4. The figure shows contours of the mass density of the 
phantom dark matter of the benchmark model seen edge-on (x-z 
plane, j/ = 0). The colour bar shows the logarithm of the PDM 
density in Mq kpc~ . In addition to the oblate phantom dark 
matter (PDM) halo, the horizontal overdensity shows clearly the 
PDM in the host disk and the vertical over density shows the 
PDM in the polar ring. The dashed circles are to demonstrate 
that the phantom halo is nearly spherical and slightly oblate. 



ters of these rings are such that the density in the center is 
zero and positive everywhere else. Note that these parame- 
ters are adopted exactly as per Combes & Arnaboldi (1996), 
but that in reality, some freedom on the mass of the stel- 
lar components in both the host and ring of NGC 4650A is 
possible. As we do not intend here to make a full detailed 
fit of the rotation curves of NGC 4650A, which will be the 
topic of a following paper, including other individual polar 
ring systems observed in HI with the WSRT (Westerbork 
Synthesis Radio Telescope), we keep the benchmark model 
as such. From this density of baryonic matter, we compute 
the corresponding PDM density using Eq. 4. The computed 
distribution of PDM in the plane orthogonal to both the disk 
and ring is plotted in Fig. 4. This figure illustrates that, in 
addition to the oblate PDM halo (isothermal at large radii), 
there are also two PDM disks aligned with the baryonic disks 
of the host and of the polar ring. 

To investigate whether the results we obtain (see 
Sect. 6) for this benchmark model are actually a generic 
prediction of Milgromian dynamics, we will then vary the 
parameters of this benchmark model in 5 different ways 
(including changing the ring into a Miyamoto-Nagai disk) 
computing a total of 45 models spanning a wide range of 
parameters. All models and their different parameters are 
summarized in Table 1. 

• First, the density of the polar rings and accord- 
ingly their mass, MpR, relative to the mass of the host 
galaxy is varied. The resulting models are collected into 
two sequences: Sequence 1 and Sequence 2. Starting from 
MpB. = 0, the ring mass is increased to Mpn. — 
(0.1, 0.25, 0.33, 0.5, 0.75, 1 and 1.45) x Mdisk. Sequence 1 has 
a constant disk mass of Mdisk = 11 x 10® M©, Sequence 2 
hasMdisk = 33x 10® Mq. 

• To obtain Sequence 3, Sequence 1 is repeated while re- 
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Figure 5. Rotation curves of the host galaxy (black, solid line) 
and polar ring (black, dashed line) for the benchmark model. 
These rotation curves are derived from the Milgromian potential. 
For comparison, the blue lower lines show the rotation curves 
derived from the Newtonian potential of the same model. Obser- 
vationally, the rotation velocity in the host is generally obtained 
indirectly from the measured stellar velocity dispersion, which 
means that the maximum velocity is most likely not measured in 
the very fiat part. To account for this issue, the theoretical rota- 
tion velocity from the models are computed at both r = 40 kpc 
and r = 15 kpc. The circle and cross at these two radii thus mark 
the rotation velocities that are summarized in Fig. 8. In the case 
where both the host and PR are gas-rich, the PR curve should 
start where the host curve ends (e.g., at 15 kpc) to avoid coUi- 
sional orbit crossing. 



placing the polar ring by a polar disk of mass Mpo- This 
polar disk is shaped like a Miyamoto-Nagai density distri- 
bution with hr = 0.748 kpc and h^, = 0.3 kpc. The polar disk 
mass, AfpD, is varied analogously to Sequence 1 and 2. 
In this sequence, the model with Mpo/Mdisk = 1 is sym- 
metric, the rotation curves in the galactic plane and in the 
polar plane are consequently identical. 

• To investigate the influence of the polar ring shape, we 
vary the shape of the ring in Sequence 4. The size parameters 
of the gaseous ring {hf^ and hf^) and of the stellar ring (/ij' 
and hf^) are summarized in Table 1. 

• Models of Sequence 5 have again the same structural 
parameters as the benchmark model, but their densities are 
scaled such that their total masses (of the whole system) 
range from 6.8 x 10^ Mq to 10^^ M0. The size parameters 
remain unchanged. 

• Sequence 6 is a series of four different Milgromian po- 
tentials computed for the benchmark model for n = 1, 2, 3, 4 
in Eq. 2, to check whether the qualitative results are inde- 
pendent of the i^-function. 



6 RESULTS 

The rotation curves of the benchmark model are presented 
in Fig. 5, both in Newtonian without dark matter and Mil- 
gromian dynamics. In both cases, the velocity is higher in 
the ring, but of course, Milgromian dynamics is needed to 
get the amplitude and shape of the rotation curves right. 
This is a key result for Milgromian dynamics, since obser- 
vational studies (lodice et al. 2003) have measured rotation 
velocities in extended polar rings being systematically larger 



than in the host galaxies when both systems are seen roughly 
edge-on. If this result is generic, it means that this particular 
property of polar ring galaxies does not exclude Milgromian 
dynamics. To investigate whether this is indeed a generic 
prediction of Milgromian dynamics, we vary the parameters 
of this benchmark model as explained in Sect. 5. Altogether, 
the rotation velocities, corresponding to the line-of-sight ve- 
locities that would be measured when observing the PRG 
edge-on, of 45 models in total are evaluated. The rotation 
curves of the benchmark model are presented in Fig. 5. Fig. 6 
contains all evaluated rotation curves of Sequences 1-4. See 
the caption for more details. Sequence 6 is plotted in Fig. 7. 
Finally, the shape of the rotation curves of Seq. 5 are all 
similar to the benchmark model, but their amplitude varies 
as in Fig. 8. 



6.1 Newtonian dynamics 

Fig. 5 and the first column of Fig. 6 show the rotation curves 
computed using standard Newtonian dynamics without dark 
matter for Sequence 1. As we can see, the velocities at radii 
larger than approximately 6 kpc are larger in the polar plane. 
The reason for this is that the closed orbits are much more 
eccentric in the host galaxy than in the plane of the polar 
ring due to the compact host galaxy and the extended ring 
(see Fig. 2), and when observing the polar ring galaxy edge- 
on, the minimum velocities of these eccentric orbits (point A 
in Fig. 2) are measured. The eccentricity can be explained 
as follows. The potential due to the compact host galaxy 
component appears nearly spherical at large^ radii for test 
particles in both the plane of the host galaxy and the plane 
of the polar ring. The potential generated by the extended 
polar ring, however, does appear spherical to particles orbit- 
ing within the ring, but not to particles orbiting in the host 
plane. ^ This gives rise to lower line-of-sight velocities when 
the host disk and polar rings are seen edge-on. 

At radii smaller than the size of the hole (the region 
where the density increases with radius), the eccentricity 
argument turns around. The test particles are near the cen- 
ter of the galaxy and accordingly near the center of the hole 
of the polar ring. At these radii, all test particles experi- 
ence the gravitational field caused by the polar ring rather 
spherically, while the potential by the galactic disk appears 
spherical only to the test objects orbiting within the disk, 
not to those moving in the polar plane. The transition ap- 
pears around approximately 6 kpc, i.e. between the radial 
size of the galactic disk, which is smaller than 6 kpc, and 
that of the polar ring, which is larger than 6 kpc. 

This ellipticity of the orbits in the host is enhanced by 
the actual presence of the hole at the center of the polar 
ring, because, in order to have the same mass in the polar 
structure as in a corresponding disk, one needs to increase 
the density at large r in the polar plane, making it more 
extended (see Fig. 3). 



* In the context of investigated PRG models, large means larger 
than the size of the host galaxy. 

^ We know this because in Newtonian dynamics the linearity of 
Poisson's equation allows us to separate and linearly add the dif- 
ferent potentials. 
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6.2 Milgromian dynamics 

In Milgromian dynamics, all investigated models (see Fig. 6) 
that (i) feature a polar ring (i.e. Sequences 1, 2, 4, 5 and 
6), tiiat (ii) has a total mass (i.e. gaseous plus stellar mass) 
comparable to the mass of the host galaxy (within a factor of 
~2), show higher velocities in the polar plane at radii larger 
than approximately 6kpc. Indeed, the reason is the same as 
in Newtonian dynamics without dark matter, and is even 
boosted by the additional gravity provided by Milgromian 
dynamics. Because the host is more compact than the ring, 
it appears more spherical to particles orbiting in the ring at 
large radii than the ring appears to particle orbiting in the 
host at the same radii. Hence, also in Milgromian dynamics, 
the closed orbits are in the galactic disk more eccentric than 
in the plane of the polar ring. This gives rise to lower line 
of sight velocities in the host for typical observed systems 
where the host and polar rings are seen approximately edge- 
on. Because Milgromian dynamics adds a disk component of 
PDM to the host and to the ring (see Fig. 4), this effect is 
even amplified compared to Newtonian dynamics. 

In Sequence 1 (Fig. 6, 1st column), one can see that de- 
creasing the polar ring mass compared to the host gradually 
cancels the above effect, because the gravity generated by 
the polar ring becomes more and more negligible when de- 
creasing its mass. For the benchmark model, at a radius of 
less than approximately 6 kpc, test stars orbiting in the po- 
lar ring have smaller velocities than those at the same radius 
in the host disk, because they are close to or in the polar 
hole and the mass enclosed by their orbits is comparably 
small. This transition radius changes to larger radii when 
decreasing the polar ring mass, to gradually arrive at the 
situation of no polar ring, where polar orbits have velocities 
systematically lower than those in the disk (equal at large 
radii) . 

The same effect is observed in Sequence 2 (Fig. 6, second 
column) for more massive systems. Note that the velocities 
in the host decrease with declining polar ring mass, due to 
the decreasing gravity of the polar ring, but less so than the 
velocities in the ring. The reason for this slower decrease is 
that the effect of decreased gravity is compensated by the 
effect of decreasing eccentricity for particles orbiting in the 
host. 

On the other hand, in Sequence 3 (Fig. 6, 3rd column), 
the rotation velocities are, at radii larger than 15 kpc, very 
similar in both planes, because of the special symmetry of 
these models. This emphasizes the role played by the hole 
at the center of the ring in the other Sequences. 

Varying the form and size of the stellar and gaseous 
holes in the polar ring in Sequence 4 (Fig. 6, 4th column) 
shows that for a host and ring of comparable mass, the effect 
is quite generic in the presence of a hole. 

The shape of the rotation curves of Sequence 5 are all 
similar to the benchmark model, but their amplitude varies 
and is discussed in Sect. 6.3. Sequence 6 is a series of four dif- 
ferent Milgromian potentials computed for the benchmark 
model for n = 1, 2, 3, 4 in Eq. 2, and confirms that the quali- 
tative results are independent of the i^-function as illustrated 
on Fig. 7. 

Let us finally note that if the host galaxy, here repre- 
sented by a Miyamoto-Nagai disk, is replaced by a disk that 
falls off exponentially and therefore much faster, the host 
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Figure 7. Rotation curves of the benchmark model in Milgro- 
mian dynamics applying different i^-functions, u„{x) with n = 
1,2,3,4 in Eq. 2. 



galaxy would appear more compact and the described effect 
would therefore be even stronger at small and intermediate 
radii, and unchanged at large radii. 



6.3 Tully-Fisher relation 

We compare our theoretically obtained rotation curves to 
the observed PRGs of lodice et al. (2003), showing, in Fig. 8, 
the maximal rotational velocities of both the host and polar 
ring. This gives rise to a luminous Tully-Fisher relation. For 
the baryonic Tully-Fisher relation (BTFR), we note that 
Milgromian dynamics predicts V^ — aoGMh for spherical 
systems (McGaugh 2011) where V is the asymptotic cir- 
cular velocity, and Mb the baryonic mass. However, po- 
lar ring galaxies are not only non-spherical but also non- 
axisymmetric objects. Because of this, there is no iron-clad 
prediction for the BTFR in such objects: our models indi- 
cate that the hosts typically exhibit asymptotic velocities 
below this value of V, while the polar rings are closer to the 
prediction, but can typically also exceed this velocity. 

In order to explore the relevant mass range in our mod- 
els, the corresponding data points in absolute B-band mag- 
nitude vs. rotation velocity are compared to the models of 
Sequence 5 (which varies the total mass), by assuming a 
mass to light ratio of M/Lb = 4 as in Combes & Arnaboldi 
(1996). We assume that the line width A«2o = W20 = 2i'max 
equals twice the maximum line of sight velocity® . 

From observations, the polar rings are generally very 
extended and feature large quantities of HI gas. The ro- 
tation curve in this polar plane can therefore observation- 
ally be measured even at large radii where it is very flat. 
Compared to the ring, the host galaxy is usually rather 
small and has relatively little gas. The rotation velocity in 
the host is generally obtained indirectly from the measured 



^ Note that the observed line widths Ai)20 are converted into 
velocities by assuming Ai)20 = 2t;niax (e.g. Verheijen 2001). De- 
pending on line width broadening effects, the actual velocities 
may be systematically smaller than Av2o/2. In the context of the 
theoretical data this would imply that NGC 4650A had actually 
a smaller mass to light ratio. 
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Figure 6. The figure presents the rotation curves of all models of Sequence 1 to 4. Solid lines refer to Milgroniian rotation velocities 
in the host galaxy, long-dashed lines to the polar plane. In Seq. 1 also the Newtonian rotation curves are shown (dotted lines: host 
galaxy; short-dashed: polar ring; no dark matter halo). Sequence 1 features a host disk, a bulge and a polar ring. While the disk mass 
A^disk = 11 X 10^ Mq is constant, the mass of the polar ring is varied from Mpa = . . . 1 .45 Afjisij (the bottom left panel shows 
the benchmark model, see also Fig. 5). Sequence 2 is similar to Sequence 1 but is three times as dense and consequently massive 
(/^disk = 33 X 10^ Mq). Sequence 3 features a host disk, a bulge and a polar disk instead of a ring. Models of this sequence have a disk 
mass of M(jig|j = 11 X 10^ Mq, the mass of their polar component is again variable. Sequence 4 features a host disk, a bulge and a polar 
ring. Models of this sequence have a fixed total mass M = 27.1 X 10^ Mq and the size parameters of their polar ring components are 
variable. 



stellar velocity dispersion (lodice et al. 2003, 2006; lodice 
2010), which means that the maximum velocity is most likely 
not measured in the very flat part and does therefore not 
equal the maximum velocity of the theoretical potential de- 
rived from the observed density distribution. To account for 
this issue, the theoretical rotation velocity from the models 
are computed at both r = 40 kpc and r = 15 kpc, where 
v{r = 15 kpc) ^ v(r — 40 kpc) (see Fig. 5). 

In Fig. 8, the circles show the rotation velocity mea- 



sured in the hosts, the arrow heads the ones measured in 
the polar rings, both for some observed PRGs and for the 
models of Sequence 5 (circles and squares as per Fig. 5). 
We see that our models reproduce fairly well the observa- 
tions, with velocities systematically larger in the rings, but 
also with comparable offsets. Since the benchmark model on 
which the Sequences are based was inspired from the pro- 
totypical PRG galaxy NGC4650A, it comes as no surprise 
that this PRG is best fitted by these models. 
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Figure 8. Comparison of observational data of PRGs witli our 
numerical results using the luminous TuUy-Fisher relation. The 
plot shows the the absolute B-band luminosity in magnitudes vs. 
the rotation velocity. Each arrow refers to one galaxy or galaxy 
model. The blue and red data points are adopted from lodice et al. 
(2003) and represent measurements of various PRGs. The circles 
show the rotation velocity measured in the hosts, the arrow heads 
the ones measured in the polar rings. The blue data correspond to 
NGC 4650A. For the theoretical data points (black), the squares 
show the rotation velocity in the host galaxy at r = 40kpc, the 
circles show the rotation velocity in the host at r = 15 kpc (where 
it is actually measured), the arrow heads point to the polar ring 
rotation velocity at r = 40 kpc. The theoretical data are obtained 
from models of Sequence 5. The absolute B-band magnitude is 
calculated from the total mass using a mass to light ratio of 
M/ L-g, = 4 like was assumed by Combes & Arnaboldi (1996). 



7 CONCLUSIONS 

The conclusion from all the investigated models, and the 
bottom line of this study, is that Milgromian dynamics nat- 
urally predicts that rotation velocities would be higher in 
the polar rings than in the hosts. This generically happens 
when the ring is more extended than the host and of com- 
parable mass and both are observed approximately edge-on. 
It does not apply to faint polar rings, or to polar rings of 
similar radial size as the host. Given the wide range of model 
parameters covered within this study, this general result ap- 
pears quite robust in Milgromian dynamics. What is more, 
the magnitude of the velocity offset predicted by the mod- 
els is also comparable to the observed one (see Fig. 8). We 
however did not attempt to precisely fit the full rotation 
curves of individual polar ring galaxies, which will be the 
subject of future work, based on observations performed on 
a sample of 10 such systems at the Westerbork Synthesis 
Radio Telescope (WSRT)'^. These and other upcoming pre- 
cise measurements of rotation curves of individual polar ring 
galaxies (e.g. the future measurements announced in lodice 
2010) should thus allow more stringent tests, and these will 
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become benchmark objects with which to test gravity in the 
coming years. 
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Scq. 



Md [10« M0] 



MpR/Md 



fcgf'^ [kpc] feff [kpc] fef^ [kpc] fef^ [kpc] 



11.0 



0.1, 0.25, 0.33, 0.5, 0.75, 1, 1.45 



6.8 



5.95 



15.3 



3.4 



33.0 



0.1, 0.25, 0.33, 0.5, 0.75, 1, 1.45 



6.8 



5.95 



15.3 



3.4 



11.0 



0.1, 0.25, 0.33, 0.5, 0.75, 1, 1.45 



0.748 



0.748 



11.0 



1.45 



15.3 


3.4 


15.3 


3.4 


7.8 


4.95 


15.3 


3.4 


5.8 


4.95 


15.3 


3.4 


6.8 


5.95 


15.3 


3.4 


6.8 


5.95 


11.3 


7.4 


6.8 


5.95 


6.8 


5.95 


10.8 


9.95 


10.8 


9.95 



2.8, 3.7, 5.6, 7.5, 8.4, 11.2, 13.0, 

14.6, 16.4, 18.4, 20.7, 23.7, 30.6, 

33.6, 36.8, 41.3 



1.45 



6.8 



5.95 



15.3 



3.4 



Table 1. Model parameters. In the presented models of PRGs, each ring is built by the difference of two Miyamoto-Nagai disks (see the 
text for more details), whose radial parameters are hn and /ir2. Model 7 of Seq. 1 and Model 4 of Seq. 4 equal the benchmark model. 



The models of Seq. 3 do not feature a polar ring but a polar disk of stars and gas with radial parameters hf-^ 
in this table, it corresponds to the benchmark model with 4 different i^-functions. 



gas/st 



Scq. 6 is not included 



Plummer H. C, 1911, MNRAS, 71, 460 
Reshetnikov V. R, Combes F., 1994, A&A, 291, 57 
Sackett R D., Rix H.-W., Jarvis B. J., Freeman K. C, 1994, 

ApJ, 436, 629 
Sackett R D., Sparke L. S., 1990, ApJ, 361, 408 
Schweizer F., Whitmore B. C, Rubin V. C, 1983, AJ, 88, 

909 
Strigari L. E., , 2012, arXiv: 121 1.7090 
Tiret O., Combes F., 2007, A&A, 464, 517 
Verheijen M. A. W., 2001, ApJ, 563, 694 
Whitmore B. C, Lucas R. A., McElroy D. B., Steiman- 

Cameron T. Y., Sackett R D., Oiling R. R, 1990, AJ, 

100, 1489 



